Quasiparticle States at a <i-Wave Vortex Core in High-T c Superconductors: Induction 

of Local Spin Density Wave Order 
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The local density of states (LDOS) at one of the vortex lattice cores in a high Tc superconductor 
is studied by using a self-consistent mean field theory including interactions for both antiferromag- 
netism (AF) and d-wave superconductivity (DSC). The parameters are chosen in such a way that 
in an optimally doped sample the AF order is completely suppressed while DSC prevails. In the 
mixed state, we show that the local AF-like SDW order appears near the vortex core and acts as 
an effective local magnetic field on the quasiparticles. As a result, the LDOS at the core exhibits 
a double-peak structure near the Fermi level that is in good agreement with the STM observations 
on YBCO and BSCCO. The presence of local AF order near the votex core is also consistent with 
the recent neutron scattering experiment on LSCO. 



PACS numbers: 74.25Jb, 74.50. +r, 74.60. Ec, 74.20.-z 



The quasiparticle states at the vortex core in the mixed 
state of a superconductor have been one of the major 
interest in condensed matter physics. For an s-wave su- 
perconductor, the energy gap opened at the Fermi sur- 
face is a constant and it was predicted long time ago 
by Caroli et al. 0j that there should exist the low-lying 
bound quasiparticle states inside an s-wave vortex core. 
This prediction was later confirmed by detailed numer- 
ical computations [|]-f§] and by STM experiments on 
NbSe2 @] although a direct observation of the discrete 
levels is yet to be performed. However, for a d-wave 
pairing state as recently established in high-T c cuprates, 
the situation becomes more complex, mostly due to the 
fact that the energy gap is closed at the nodal direc- 
tion on the essentially cylindrical Fermi surface. In an 
earlier study by Wang and MacDonald based on a lat- 
tice model ||, it was shown that the local density of 
states (LDOS) at the G?-wave vortex core exhibits a sin- 
gle broad peak at zero energy. Recent low temperature 
STM experiments on YBa 2 Cu 3 7 _ 5 (YBCO) 0(a) and 
Bi 2 Sr 2 CaCu 2 8+( 5 (BSCCO) g both observed a double- 
peak structure around zero bias in the local differential 
tunneling conductance at the vortex core center, which 
identifies the widely split core states at energies ±5.5 
meV and ±7 meV, respectively. The discrepancy be- 
tween the theory and the experiment stimulated further 
theoretical studies J9|-[l3| on the quasiparticles in the vor- 
tex core of high-T c cuprates. Franz and Tesonovic |Tl| l 
proposed an explanation of the observed double-peak 
structure in terms of a mixed d x 2 +y 2 +id xy pairing state. 
Such a pairing state was previously suggested Q to be 
realized through the field-induced second phase transi- 
tion as motivated by the observation of a plateau in ther- 
mal conductivit y jig ]. The origin of this plateau is still 
hotly debated ]16|. By pointing out that the parame- 



ter values are unrealistically chosen in Ref . 1 1 that the 



d x 2 +y 2 + id xy state even already exists in zero magnetic 
field, Yasui and Kita applied the "Landau-level ex- 



pansion method" to study the quasiparticles in d-wave 
vortex lattice based on the continuum model. They found 
that the double-peak structure may be considered as in- 
herent to systems with short coherence length. Since 
the splitting sensitively depends on the field strength, 
the validity of this scenario still needs to be clarified 
theoretically and experimentally. On the other hand, 
the superconducting vortex with antiferromagnetic (AF) 
core was also predicted based on either the SO(5) the- 
ory JlTj or the standard t-J model with spin-charge sep- 
aration |Q, which leads to a featureless LDOS and is in 
disagreement with the STM data on the optimally doped 
YBCO §(a) and BSCCO @. Partly motivated by the 
observation of magnetic vortex cores in recent neutron 
scattering experiment by Lake et al. [l9| ] on optimally 
doped La 2 _ a; Sr a ;Cu04, we present in this Letter an alter- 
native mechanism for the double-peak structure in the 
LDOS around zero energy. We show that due to electron 
correlations, the AF-like spin density wave (SDW) order 
can develop locally around the vortex core and vanish 
in the superconducting regions. The lift of the spin de- 
generacy leads to the splitting of the zero energy peak. 
The induced SDW order around the vortex core man- 
ifests the repulsive electron interaction responsible for 
the strong spin fluctuations in the underdoped region 
of high-r c cuprates. In fact, the coexistence of super- 
conducting (SC) and SDW orders has been theoretically 
studied |2fj|-p3||, which shows a rich phase diagram with 
a classic AF order at half filling, striped phase at under- 
doping, a d-wave SC at optimal doping. 

We start with a generalized Hubbard model defined on 
two dimensional (2D) lattice. By assuming that the on- 
site repulsion is solely responsible for the antiferromag- 
netism while the nearest neighbor attraction causes the 
d-wave superconductivity, we can construct an effective 
mean-field (MF) model |2j| to study the vortex physics 
in the mixed state: 
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Here c- la annihilates an electron of spin a at site i. 
The summation is over the nearest neighbor sites, /i is 
the chemical potential. m\, a = U{a Ci a ) is the spin- 
dependent Hartree-Fock potential at site i, where U is 
the strength of on-site repulsion. Ay = y (c^Cji — Cij.Cj-[-) 
is the spin-singlet d-wave pair potential, where V is the 
strength of nearest neighbor effective electron-electron 
attraction. As a phenomenological model, we do not 
intend to address the microscopic mechanism for this 
attraction. In the mixed state, the magnetic field ef- 
fect was included through the Peierls phase factor tpy — 
|^ f*' A(r) • dr, where $o = hc/2e is the superconduct- 
ing flux quantum. By assuming the superconductor un- 
der consideration is in the extreme type-II limit where 
the Ginzburg-Landau parameter k = A/£ goes to infin- 
ity so that the screening effect from the supercurrent is 
negligible. Therefore, the vector potential A can be ap- 
proximated by the solution V x A = Hz where H is 
the magnetic field externally applied along the c axis. 
The enclosed flux density within each plaquette is given 
by So Vij = 2lT $ Q a ■ A similar mean-field Hamiltonian 
can also be arrived at within a t-U-J model proposed 
recently |24| . 

We diagonalize the Hamiltonian Eq. (|l|) by solving the 
BdG equation: 

?0fr -&)(*)-*(*)■ • 

where (u^, vf^) is the quasiparticle wavefunction corre- 
sponding to the eigenvalue E n , the single particle Hamil- 
tonian 7iij,(j = — te zip ' 5 ^i + s j + ( m i.s ~ M^y- Notice that 
the quasiparticle energy is measured with respect to the 
Fermi energy. The self-consistent conditions read: 

m ia = Uj2\<\ 2 f(En), (3) 

n 

and 




where the Fermi distribution function f(E) = 
l/(e E / klsT + 1). Here the summation is also over 
those eigenstates with negative eigenvalues thanks to 
the symmetry property of the BdG equation: If 
(u§,v£ l ,v§,v? l ) Tranapoae is the eigenfunction of the 4 x 
4 equation in the spin space with energy E n , then 
(v§*,-v? l *,u§*,-u§*) Trans P°° e up to a global phase 
factor is the eigenfunction with energy —E n . 



Hereafter we measure the length in units of the lat- 
tice constant a and the energy in units of the hopping 
integral t. Within the Landau gauge the vector poten- 
tial can be written as A = (—Hy, 0, 0) where y is the 
y-component of the position vector r. We introduce the 
magnetic translation operator T mn r = r + R where the 
translation vector R = mN x e x + nN y e y with N x and 
N y the linear dimension of the unit cell of the vortex 
lattice. To ensure different T mn to be commutable with 
each other, we have to take the strength of magnetic 
field so that the flux enclosed by each unit cell has a 
single-particle flux quantum, i.e, 2$ - Therefore, the 
translation property of the superconducting order param- 
eter is A(7^ rm r) = e JX ( r ' R * 1 A(r) where the phase accu- 
mulated by the order parameter upon the translation is 
x(r,R) = |~A(R)-r— 4mmr. From this property, we can 
obtain the magnetic Bloch theorem for the wavefunction 
of the BdG equations: 

( U k . a (%nnr) \ _ ik-R ( e^^^U^r) \ 

Here f is the position vector defined within a given 
unit cell and k = m^n^ &x + ^ N &v witl1 m ^y 

( ). I M x , y — 1 are the wavevectors defined in the first 

Brillouin zone of the vortex lattice and M X N X and M y N y 
arc the linear dimension of the whole system. The vortex 
carrying the flux quantum hc/2e is the generic feature of 
the pairing theory for superconductivity. Therefore, it is 
not surprising that in the slave boson approach to the 
t-J model, the vortex always carries hc/2e flux quantum 
if the magnetic field is assumed to act on electrons 
only through the spinon degrees of freedom. 

As a model calculation, we take the following pa- 
rameter values: The pairing interaction is V = 1.0, 
and the filling factor, which is defined as nj = 
Si a{ c \<j c itr) /N x N y with the summation over one unit 
cell, is fixed to be 0.84 so that the chemical potential 
needs to be adjusted each time the on-site repulsion 
U is varied. For our interest in the low energy quasi- 
particle states, we only consider the zero temperature 
limit. We have typically considered the unit cell of size 
N x x N y = 42 x 21, and the number of the unit cells 
M x x My = 21 x 42. This choice will give us a square 
vortex lattice. We use exact diagonalization method to 
solve the BdG equation (|2|) self-consistently: To allow 
the inhomogeneity of all physical quantities, randomly 
distributed Ay and m\u are taken as initial parameters; 
the newly obtained Ay and m\^ a are then substituted 
back into the equation. The above procedure is repeated 
until the convergence with required accuacy is achieved. 
In the absence of magnetic field, we have reproduced 
the results reported in previous work |23f| including an 
AF SDW order, a stripe phase, and a d-wave SC phase 
when the system is doped away from the undoped to 
the optimally doped region. In the present work, we are 
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mainly concerned with the electronic structure around 
the vortex core in the optimally doped region. In this 
region, the SDW order is strongly suppressed and the d- 
wave SC order is homogeneous in real space. However, 
when a magnetic field is applied to drive the system into 
the mixed state so that the <i-wave order parameter is 
suppressed around the vortex core, we find that as the 
on-site repulsion is increased to about 1.5, the SDW or- 
der is nucleated around the vortex core. Typical results 
on the nature of the vortex core is displayed in Fig. |l| 
with the on-site repulsion U = 2. As shown in Fig. |l](a), 
each unit cell accommodates two superconducting vor- 
tices each earring a flux quantum hc/2e. The d-wave SC 
order parameter vanishes at the vortex core center and 
starts to increase at the scale of the coherence length 
£o to its bulk value which is about 0.1 for the chosen 
parameter values. Fig. 0(b) displays the spatial distri- 
bution of the staggered magnetization of the local SDW 
order as defined by M s = (— 1)'5* with SI — n-^ — n^. 
Clearly, the maxima strength of M s appears at the vor- 
tex core center and decays also with a scale of £o to 
zero into the superconducting region. More interestingly, 
the SDW order parameters has opposite polarity around 
two nearest neighbor vortices along the x direction. We 
have compared the free energy between this configura- 
tion with that obtained by switching the orientation of 
the SDW order around one of the two nearest neighbor 
vortices, and found that the free energy for both config- 
urations is very close but it is always lower in the former 
case. Therefore, the induction of the SDW order around 
the vortices reduces the four-fold rotational symmetry 
of the whole system to the two-fold, and the period of 
the translational symmetry of the vortex lattice along 
the x is doubled. This result is understandable when we 
notice the zero-field result : The homogeneous super- 
conducting order in the optimally doped region derives 
from the melting of the strongly overlapped quasi-one 
dimensional (Q1D) superconducting stripes (i.e, soliton- 
like AF anti-phase domain boundaries, at which the AF 
SDW order changes sign). The development of the Q1D 
stripes breaks at the beginning the four-fold rotational 
symmetry of the system. Therefore, it seems that the 
development of the local SDW order around the vortices 
in the optimally doped region can be regarded as a du- 
ality of the development of the local SC order around 
the AF stripe in the underdoped region. On the other 
hand, the appearance of the SDW order around the vor- 
tex strongly affects the electron density n\ = J2a ni ^- 
As shown in Fig. |l|(c), at the vortex core center, where 
the SDW amplitude reaches the maximum, the electron 
density is strongly enhanced and is very close to unity, 
which is characteristic of the bulk AF-like SDW order 
at the half filling. Therefore, the hole charge density is 
depleted in the vortex core center. The depletion of the 
hole charges near the vortex core center is compensated 
by the corresponding enhancement along ir/4 and 37r/4 



directions with respect to the underlying crystal lattice, 
which correspond to the nodal directions on the Fermi 
surface. This kind of charge inhomogeneity is closely re- 
lated to the development of the local SDW order around 
the vortex core due to the large on-site repulsion. Finally, 
with the chosen parameter values, we also find that when 
the on-site repulsion U is increased to 3.5, the SC vortex 
state is collapsed and the SDW order becomes dominant. 
Our numerical analysis seems to be consistent with the 
recent argument p5| of the magnetic field driven quan- 
tum phase transition from the SC state into a state with 
microscopic coexistence of SC and SDW orders to un- 
derstand the recent neutron scattering experiments [fl9| . 
Notice that the quantum fluctuation of the SDW order 
parameter, which plays an important role in addressing 
correctly the spin excitations, has been neglected in the 
present model because our interest is in the quasiparticle 
states. 

We now turn our attention to the quasiparticle state 
at the vortex core center. The local density of states is 
defined by 



Pi(E) 
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(6) 



where f'(E) is the derivative of the Fermi distribution 
function. pi(E) is proportional to the local differential 
tunneling conductance which could be measured by STM 
experiments 26 1. In Fig. || we plot the LDOS as a func- 
tion of energy at the vortex core center for different values 
of on-site repulsion. For comparison, we have also dis- 
played the LDOS at the midpoint between two nearest- 
neighbor vortices along the x direction, which resembles 
that for the bulk system. The asymmetry line shape 
in pi(E) with respect to zero energy reflects the lack of 
particle-hole symmetry as the chemical potential p de- 
viates from zero for n/ being less than the half filling 
(rif = 1). As can be seen from Fig. §(a), when U = 
for which no local SDW order is induced, the LDOS at 
the core center shows a single resonant peak around the 
Fermi energy, which is similar to that reported by other 
authors ||. When U is sufficiently large that the local 
SDW order develops around the vortex core, the LDOS 
peak at zero energy is split into a double-peak struc- 
ture (see Fig. §(b)). The splitting comes from the fact: 
When the SDW order is localized around the vortex cen- 
ter, the spin-dependent potential, which can be rewrit- 
ten as TOi,|(x) = U(n,i ± SD/2, plays the role of a local 
magnetic field interacting with the electrons via the Zee- 
man coupling. When U is increased, S* at the vortex 
core center is enhanced, and the combination of them 
enlarges the Zeeman interaction. As a consequence, the 
LDOS peak at zero energy is further split (see Fig. |^(c)). 
This splitting of the LDOS at the vortex core center is in 
good agreement with the STM experiments on the opti- 
mally doped YBCO and BSCCO. The induction of the 
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local AF order at the vortex core is consistent with recent 
neutron scattering experiment on LSCO Jl9| . 
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FIG. 1. The amplitude distribution of the d-wave SC order 
parameter |Ad| (a), the staggered magnetization M s (b), and 
the electron density rii (c) in one magnetic unit cell obtained 
at the zero temperature. The size of the cell is 42 x 21. The 
strength of the on-site repulsion U = 2. The other parame- 
ter values: The d-wave pairing interaction V = 1; the filling 
factor rif = 0.84. 

FIG. 2. The zero-temperature LDOS at the vortex core 
center (red-solid line) with various strength of on-site repul- 
sion U = (a), U = 2 (b), and U = 3 (c). Also displayed is the 
zero-temperature LDOS at the midpoint (green-dot-dashed 
line) between two nearest neighbor vortices along the x direc- 
tion. The other parameter values are the same as Fig. hi 



C. Caroli, P. G. de Gennes, and J. Matricon, Phys. Lett. 
9, 307 (1964). 

F. Gygi and M. Schluter, Phys. Rev. Lett. 65, 1820 
(1990); Phys. Rev. B 41, 822 (1990); 43, 7609 (1991). 
J. D. Shore et ai, Phys. Rev. Lett. 62, 3089 (1989). 
Yu-Dong Zhu, F.C. Zhang, and M. Sigrist, Phys. Rev. B 
51, 1105 (1995). 

H. F. Hess et al, Phys. Rev. Lett. 62, 214 (1989). 

Y. Wang and A.H. MacDonald, Phys. Rev. B 52, R3876 

(1995). 

(a) I. Maggio-Aprile et al, Phys. Rev. Lett. 75, 2754 
(1995); (b) Ch. Renner et al, Phys. Rev. Lett. 80, 3606 
(1998). 

S. H. Pan et at, Phys. Rev. Lett. 85, 1536 (2000). 
A. Himeda et al, J. Phys. Soc. Jpn. 66, 3367 (1997). 
Y. Morita, M. Kohmoto, and K. Maki, Phys. Rev. Lett. 
78, 4841 (1997); 79, 4514 (1997); M. Franz and M. 
Ichioka, ibid. 79, 4513 (1997). 

M. Franz and Z. Tesanovic, Phys. Rev. Lett. 80, 4763 
(1998). 

M. Takigawa, M. Ichioka, and K. Machida, Phys. Rev. 
Lett. 83, 3057 (1999). 

K. Yasui and T. Kita, Phys. Rev. Lett. 83, 4168 (1999). 

R. B. Laughlin, Phys. Rev. Lett. 80, 5188 (1998). 

K. Krishana et al, Science 277, 83 (1997). 

H. Aubin et ai, Phys. Rev. Lett. 82, 624 (1999). 

D. P. Arovas et ai, Phys. Rev. Lett. 79, 2871 (1997). 
Jung Hoon Han and D.-H. Lee, Phys. Rev. Lett. 85, 1100 
(2000). 



B. Lake et ai, Science 291, 1759 (2001); cond- 



mat/0104026 



L. Balents, M. P. A. Fisher, and C. Nayak, Int. J. Mod. 
Phys. B 12, 1033 (1998). 

S. Sachdev, C. Buragohain, and M. Vojta, Science 286, 
2479 (1999). 



J. H. Han, Q.-H. Wang, and D.-H. Lee, cond 



mat/0006046 



I. Martin et al. , |cond-mat/0009067 



S. Daul, D. J. Scalapino, and S. R. White, Phys. Rev 
Lett. 84, 4188 (2000). 



E. Demler, S. Sachdev, and Y. Zhang, cond 



mat/0103192 



M. Tinkham, Introduction to Superconductivity (McGraw 
Hill, New York, 1975). 



4 



FIG. 1(a) Zhuet al. 




FIG. 1(b) Zhuet al. 
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